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ABSTRACT 


The feasibility of performing nondestructive analysis of spent 
power reactor fuel by means of gamma ray spectroscopy was esta- 
blished. Dose rates to be expected from spent power reactor fuel 
pins offered by the AEC for study were caJculated, Parallel cal- 
culations were performed for MITR fuel elements and the results 
compared with measured dose rates of the MITR fuel elements. 
fama ray Spectra were obtained from MITR fuel elements and used 
to estimate the expected gamma ray spectra for the spent pins. 
Calculations were performed to determine the number of Pu-239 
Biesions which had occurred in the spent fuel pins and the re- 
Mme used to Calculate the expected change in the spectra due 
to the presence of Pu-239 fission product gamma rays. Recommendations 
are made to enhance the ability of M.I.T. researchers to perform 
the subject experiments on actual spent power reactor fuel. 
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i tN TRODUGELON 


The possibility of utilizing gamma ray spectroscopy for 
@etermination of £1SSlo0n product content for use in ascertain- 
mg the irradiation history of spent fuel has been studied by 
many investigators (1-4). Sovka (5) utilized the excellent 
energy resolution of Ge(Li) detectors to develop a method for 
determining spatial distribution of absolute neutron flux, 
neutron exposure, total U-235 burnup and irradiation time of 
the fuel. This method was unique in that it did not require 
absolute calibration of the detecting system since it used 
ratios of fission product activities. However, because it 
does use ratios of fission. product activities, more than one 
fission product's gamma ray(s) must be resolved. This work 
was performed using the spectra obtained from spent MITR fuel. 
Since MITR fuel has a pre-irradiation enrichment of 93%, it 
differs considerably from that of the more widely used and 
economically more important PWR or BWR power reactor fuel 
which usually has an enrichment of 2-4%. Because of the im- 
portance of accurate knowledge of the present status and past 
history of partially burned or spent power reactor fuel from 
the viewpoint of fuel cycle economics, reactor core physics, 
and most recently, safeguards, it is quite important that ad- 
vanced spectrometry methods be extended to these fuel types. 
This extension is in part, one of the aims of the M.I.T. Re- 
eeeor Physics Project which is sponsored by the AEC and under 
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whose aegis the present work was performed. For example, the 
power reactor fuel with its high percentage of U-238 would 
have a Significant number of Pu-239 fissions, and thus the 
spectra from spent power fuel would be altered because of the 
a@ifference in the fission product yields for U-235 and Pu-239, 
Self-shielding effects upon the emitted spectra would also 
obviously differ. The purpose of this study was therefore to 
determine the feasibility of performing gamma ray spectroscopy 
on spent power reactor fuel, 

In order to establish feasibility, gamma ray spectra were 
obtained from spent MITR fuel to measure the energy resolution 
of the detecting system, to obtain data from which to estimate 
the relative intensities of gamma rays present in power reactor 
fuel spectra, and to determine modifications to the scanning 
Syctem necessary for the scanning of power reactor fuel by 
Subsequent investigators. 

Calculations of the dose rate to be expected from the 
spent fuel offered by the AEC for study oe made in order to 
develop. safe handling procedures for the highly radioactive 
fuel, which was in turn necessary in order to obtain AEC 
licensing for the subject experiments. Parallel calculations 
were performed on spent MITR fuel and comparisons made with 
observed dose rates in order to check the fidelity of the cal- 
culations and to normalize the predicted results to physically 
measured quantities. Finally, calculations of the relative 


gamma ray intensities present in the power reactor fuel spec- 
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tra after a 2.75 year cooling period were made and their re- 
sults used to determine the suitability of the proposed fuel 
for analysis. Each of the above aspects are reported upon in 


detail in tne three chapters which follow. 
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Liz BOSE RATE CALCULATION FOR SPENT FUEL 


met Ilntroduction 


Prior to obtaining spent power reactor fuel it was nec- 
essary to develop handling procedures and to secure AEC li- 
cenSing approval, which necessitated estimation of the dose 
rate eecien from the spent fuel. Two fuel pins which had 
been tested in the Dresden power reactor were made available 
by the AEC for the studies. Some of the known and calculated 
information (6) on the fuel pins, numbered A4 and A28 in the 


Dresden core loading is listed in Tables 1 and 2. 
2.2 Background Information 


Many studies (7-11) have been published on the activity 
of spent fuel ranging from experimental studies based on calor- 
imetric methods to theoretical calculations. The experimental 
studies present only total activities and as such are unsatis- 
factory for dose rate claculations, since the energy depen- 
dence of thetactivity must be known, because the mass absorp- 
MicMmCcOctiteient, the bulldup factor and flux required for a 
given dose rate are all strongly dependent on energy. The 
theoretical compilations have utilized fission product data 
to calculate the activities for varying irradiation and decay 
times. Blomeke and Todd (12) have calculated the activities 
of specific fission products, and while this study thus gives 


the complete energy dependence of the spectra, the number of 
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fission products makes this method too cumbersome for estimat- 
ing the dose rate from spent fuel. Most studies divide the 
energy spectrum of fission product gamma rays into two or more 
energy groups and sum all gamma rays whose energies fall with- 
in these groups. If a sufficient number of groups is used the 
calculations will account for the energy dependence mentioned 
previously without the necessity of accounting for each fission 
product's activity on an individual basis, 

The study by Prawitz et al (7) was chosen to estimate the 
activity of the spent fuel because a sufficient number of 
energy groups are used and the decay of the activities is pre- 
sented for cooling times up to 20 years, whereas many studies 
Seenoc report the decay for times longer than one year, ‘The 
study gives the activity of fission products as a function of 
both irradiation time and cooling time and some indication is 
also given of the trend for flux variations. This latter para- 
meter is important only for fission products which have high 
neutron absorption cross sections. Statistically, this group 
of fission products represents only a small fraction of the 
Peters SLOne product activity, Li the flux dependence must 
be known the comprehensive study of Blomeke and Todd (12) 
presents the activity of spent fuel for various irradiation 
and cooling times as a function of neutron flux. While this 
study does give the dependence of activity on flux it was felt 


this advantage was offset by the fewer number of energy groups 
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used, 
m3 Method of Calculation 


In order to ensure a conservative estimate and to sim- 
plify the calculations, it was assumed the burnup occurred 
during a one year irradiation period, followed by a cooling 
time of 700 days. 

The method of calculation will be illustrated for fuel 
pin A28 with a 4 foot water shield. 

1) The 15 energy groups and their energy range are 
listed in columns 1 and 2 of Table 3. The effective energy of 
each group was assumed to be equal to the midpoint energy of 
each of the groups and this energy is listed in column 3 of 
Table 3. This energy was used for determining the values of 
the mass absorption coefficient, buildup factor and the flux 
to dose rate conversion factor to be used with the various 
energy groups. The release of an amount of energy equal to 
the effective energy was also assumed to represent one dis- 
integration. The specific energy release rates for the var- 
10us energy groups for an irradiation time of one year and 
700 days cooling time are listed in column 4 of Table 3, The 
specific energy release rates were converted to specific 
activity by dividing by the effective energy for each group 
to convert from Mev/gFP minute to disintegrations/gFP second. 
nts Specific activity is listed in column 5 of Table 3. 
Specific activity was converted to Activity by multiplying by 
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the burnup experienced by the fuel pin. 

2) Once the activities of the energy groups are known, 
Eney must be converted to fluxes present at the point in ques- 
tion. Reference (13) contains expressions relating flux to 
activity, source geometry and shield thickness. For a cylin- 
drical source of finite length the following expression is 


given for the geometry shown in Figure Il: 


\ 
o& 
S 
ye 


ae Dn) F(6, bo) (1) 


where: 


@ = undirected flux (y rays/om-sec) 
S _ = volume source strength (y rays/cm>-sec) 


Ba-eoUMLupetactol tO COnLect for multiple 


scattering 
r = radius of cylindrical source (cm) 
a= distance from source to observer (cm) 
| zZ = effective self-absorption distance (cm) 
6 -b 


F(q} bo) = fe Seqrodo 
O 


ee ez 


2 a S 
eam eer 
H, = Mass absorption coefficient for shield 
ty Sclimewneccue: Shi1eld 
a absorption coefficient of the source 


Because the diameter of the fuel pins is small the self-absorp- 
tion distance is small compared to the shield thickness, and 


ay = 
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the fuel pins may be approximated as line sources of finite 
length, 


In this case, Equation (1) simplifies to: 
ie 
27a 


Sy = line source strength (y rays/cm) 


In Table 4 the factors Hy 7b 


where: 


em (00>) 7 BENemb. ) Vand 


BF (68 b ) 
_ eas are tabulated for a 4 foot thick water shield for 
27a 


PSenmeonmeme energy groups. Values for #. and F(0,b,) were ob- 


i 
tained from Reference (13). Values from B were obtained from 
Reference (14) for a point isotropic source in water. The 
value of @ used was 40°, 

In Table 5 column 2, the values of the line source strengths 


for each group are tabulated, having been calculated from the 


relation: 


moet hey, PerCervaty (3) 
Pin Length 40 (2.54) 


The gamma fluxes received at the point 4 feet from the 
fuel pin were then calculated from Equation (2) and the values 
are contained in column 3 of Table 5. The flux necessary for 
1 mr/hr obtained from Reference (15) is tabulated in column 4 


Of Table 5 for each energy group. The dose in mr/hr for each 
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MOU tomedouEdeoa tm Column 5 OF Table 5. The total dose is 
mien sowmnron the doses from each energy group. 
Cateutatienouy-reosalso performed for a 2 foot thick water 
shield and in 1 meter of air. The dose rates for fuel pin A4 
were calculated from those of A28 by MU Lerplly Ine yenenee ra cle 
of the burnup received. The results of these calculations are 
Summarized in Table 6. 
iimonucem sto cmeck the validity of the calculations a com- 
Pattson was made between values of dose rate calculated for 2 
MITR fuel elements and observed values for these elements. 
Mne results are tabulated in Table 7. It was felt that the dis- 
crepancy was due in part to the approximation of the MITR ele- 
Meniteas a line source, Because of the relatively large size 
Of these elements appreciable attenuation would occur in the 
element. To check this assumption the calculation was perform- 
ed again for MITR eieneut 2M15 using the finite cylinder ex- 
pression, Equation (1). The radius was assumed equal to one 
half the thickness of the element, 1.5 inches. The calculat- 
ed dose rate using Equation (1) was 37,2 pao ae Conpanued= to 
43.6 mr/hr for the line source expression. Another reason 
mor the discrepancy between the calculated and observed. 
doses for the MITR elements is the intermittent operation 
of the M.I.T. Reactor, since Prawitz (7) assumed continuous 
operation in the development of the relations used in these 


calculations. 
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TABLE 6: CALCULATED DOSE RATES FROM FUEL PINS (MR/HR) 


FUEL PeMErER ZF Cees 0 ANE ts, ti. 0 
PIN AIR 


2 2 
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TABLE 7: COMPARISON OF CALCULATED AND MEASURED DOSE RATES 
BOte Vel R wr UP nL BEEMENIS, 42 -F1, H,0 SH LEED 


TOTAL2 IRRADIATION? COOLING@ CALCULATED MEASURED 
ENERGY Pin LOD PERIOD DOSE RATE DOSE RATE 
MWH (mr/hr) CubayAaaben) 


850 days 


130 days 





(a) 


Reference (2])). 


me 7 ae 





Because of the reasonably good agreement between the 
calculated and observed dose rate values for the MITR fuel it 
was felt the parallel calculations for the Dresden fuel would 
also yield adequate and even more important, conservative 
estimates of their expected dose rates. In particular, it 

“should be noted that the expected dose rates from the Dresden 
fuel pins are well below dose rates now experienced from MITR 
fuel elements under similar shielding conditions. This is an 
important conclusion because it implies that safe handling of 
the Dresden fuel can be accomplished using available MITR 
facilities and fuel handling procedures, It is therefore 
recommended that licensing eee be requested from the AEC 
on this basis. 

In view of this conclusion, it was felt that no further 
work was required on this aspect of the problem, and attention 
could therefore be focussed on the actual spectrometry, as is 


to be discussed in the following chapter. 
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Ill: EXPERIMENTAL INVESTIGATIONS 


fl ee introduction 


In order to use gamma ray spectra for analysis of spent 
fuel, certain fission product gamma rays must be resolvable 
both in intensity and in energy so that the measured count 
rate is a statistically significant measure of the activity due 
solely to the particular fission product. In order to use the 
method of analysis developed by Sovka (5) it is desirable to 
meamaote tO resolve the gamma rays of Cs-13/, Cs-134, Zr-95 and 
Pr-144, These fission products have high fission yields and 
relatively high energy gamma rays; pertinent data concerning 
these fission products is presented in Table 8. 

The experiments described below were performed on MITR 
Spent fuel in order to: 

1) Obtain data from which to estimate whether the gamma 
rays for the above fission products will be resolv- 
PpleomOuscOOlimgr times On the order of 2 1/2 years 
for the power reactor spent fuel pins described in 
Chapter 2, 

2) Determine the improvements in resolution of the spent 
fuel gamma spectra resulting from the improved re- 
Scelucionvouspresent Ge(Li) detectors, and use of the 
4096 channel analyzer, and thus ascertain if addition- 
al fission product gamma rays were suitable for analy- 
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3) Check operation of the scanning system and make alter- 
ations as necessary to enable the equipment to scan 


eideholds Spent power reactor fuel pins. 
3.2 Experimental Equipment and Procedures 


3.2.1) MITR Fuel Elements 

Mitte evemecletenes sane of the Standard MTR plate 
type construction. A detailed description is contained in 
Reference ( 5). The important feature for present purposes 
Peomienat the elements contain approximately 160 grams before 
irradiation of highly enriched uranium (93% U-235). Gamma ray 
Spectra were recorded for two elements, 7M9 and 5M22. Some 
Of the pertinent information for the elements are listed in 


maple 9. 
3.2.2) Scanning Equipment and Procedures 


Scanning equipment developed by Sovka (5) for 
scanning spent MITR fuel was used to obtain spectra from 
irradiated MT'TR fuel. A schematic diagram, shown in Figure 2, 
illustrates the major features of the system which allows the 
Scanning of spent fuel to be accomplished in the MITR spent 
fuel storage pit. The detector dewar is positioned on the 
upper movable carriage and is surrounded by lead shielding to 
reduce background. The upper movable carriage is positioned 
on the lower movable carriage, which permits scanning in two 
perpendicular directions in the horizontal plane. Extending 


from the upper carriage through the water to the fuel element 
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ny Cetin. dir rallied Eube which acts as a collimator 
and allows a well-collimated beam of gamma rays to reach the 
detector, Additional collimators made of lead with diameters 
meom 1/2 zn. to 1/8 in. may be placed in the beam of the top 
carriage to reduce the count rates to acceptable levels and 
provide additional shielding against background. [In the pre- 
sent experiment, lead sheets were also placed in the beam to 
absorb many of the low energy gamma rays in the beam while 
allowing the higher energy gamma rays to pass through rela- 
tively unattenuated. Count times of 80 to 160 minutes were 
used depending on the count rate. 

A schematic diagram of the electronic components used 
together with important settings are shown in Figure 3. A 
list of these components is shown in Table 10. The settings 
shown in Figure 3 resulted in an energy decrement of approxi- 


mately .66 keV/channel. 
3.3 Experimental Results 


Identification of the photopeaks noted in the spectra 
obtained from the spent fuel elements was made utilizing the 
energy calibration of Reference (5). The improved resolu- 
tion of the Ge(li) detectors now manufactured at M.I.T. allow- 
Sa better resolution of the spectra, while their increased 
Size and hence greater efficiency made shielding against back- 
ground more difficult. Orginally, a detector with an active 


volume of 30 com> was used, but even with 4 inches of lead 
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TABLE 10: 


EQUIPMENT 
Preamplifier 

Bias Supply 

Amplifier 

4096 Channel Analyzer 
Spectrum Stabilizer 
Memory Unit 


Readout Control 


Piece EROn TC EOULPMENT DATA 


MANUFACTURER 
Canberra Instruments 
Canberra Instruments 
Eve Waeta ap 

Nuclear Data 
Nuclear Data 
Nuclear Data 


Nuclear Data 
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‘MODEL NO. 


1408C 
3001 
1417 
ND-l161F 
ND-502 
ND-160M 


ND-160R 





shielding on the bottom of the detector and 2 inches of lead 
shielding on the side of the detector,dead times of 30-35% 

were noted even with the collimator opening also shielded by 
4 inches of lead, Additional shielding could not be added as 


the amount of shielding material was limited by the mech- 


Soe 
anical strength of the scanning equipment. This larger detector 
was replaced with a smaller detector having an 8 com> Volume 
The resulting dead times were found to be on the order of 

30% with the COllMMecorseube Open and an additional 1/8 in. 


aperture collimator inserted immediately below the detector 


when element 7M9 was placed in the holder. 
3.3.1) Fission Product Gamma Ray Spectra 
3,.3.la) Fuel Element 7M9 


The spectrum obtained from MITR fuel element 7M9 
after a cooling period of 79 days is shown in Figure 4, The 
spectrum was obtained with a 1/8 in. diameter aperture colli- 
mator for a counting time of 80 minutes. In addition three 
1/4 in. lead sheets were placed in series in the beam to pre- 
ferentially reduce low energy gamma rays in the beam. The 
spectrum shows resolved peaks for gamma rays from the decay 
of Nb-95, Pr-144, Cs-134 and Cs-137. Except for the peaks due 
to Pr-144, all resolved peaks lie between 500 and 800 keV, and 
an enlargement of this portion of the spectrum is shown in 


PIGUre 5. 
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Soe uel Element 5M2 7 


The spectrum from MITR fuel element 5M27 after a 
cooling period of 272 days is shown in Figure 6. The spec- 
trum was obtained using a 1/8 in. diameter aperture colli- 
mator for a counting time of 160 minutes. One 1/4 in. lead 
sheet was placed in the beam to preferentially reduce the low 
energy gamma rays in the beam. The spectrum shows resolved 
peaks for gamma rays due to the decay of Nb-95, Zr-95, Pr-144, 
Cs-134 and Cs-137/, The spectrum between 600 key and 800 kev 


i> Shown in Figure 7. 
oe Calculations of Gamma Ray Intensities 


Intensities of gamma rays resolved in the spectra of 
Figures 4 and 6 were calculated by finding the area, express- 
ed in terms of counts, under each peak. These calculations 
were performed using GAMANL, Reference (16), a computer code 
which identifies peaks in spectra and determines their cen- 
ters and areas. GAMANL utilizes Fourier transforms to smooth 
mene data, which eliminates much of the random Fluctdat rons si) 
Ene spectra without affecting the resolution. The background 
is then calculated and subtracted from the smoothed data and 
the resulting peaks are then analyzed for peak centers, areas 
and FWHM's, among other things. The areas are calculated in 
two ways, in the first method all the counts are added be- 


tween successive zeroes in the smoothed data with background 
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subtracted. Tne second method uses a Gaussian fit to the 
peaks and calculates the area using the height of the peak and 
its standard deviation, 

Results for the spectra shown in Figures 4 and 6 are 


shown in Tables 11 and 12 respectively. 
3.5 Conclusions 


Based upon the results obtained in the work reported in 
this chapter, the following conclusions can be drawn pertinent 
to the projected future experiments on power reactor fuel: 

1) The improved resolution of the Ge(Li) detectors now 
manufactured at M.I.T. and the 4096 channel and analyzer are 
evident in partial resolution of the Zr-95 Nb-95 doublet at 
758 and 766 keV. This consideration is not too important for 
cooling times on the order of one year or less since the peaks 
necessary for analyses are all strong peaks and tend to mask 
weaker peaks in their Compton background. After several 
years cooling, however, the peaks due to gamma rays from Pr-144 
Baltcowes-!34 will have decreased considerably due to their rela~- 
tively short half lives of 285 days and one year respectively. 
The peas due to Nb-95 and Zr-95 will decrease even more 
rapidly due to their effective half life of 65 days, and will 
be a relatively small peak after a cooling time on the order 
of 1 1/2 years. The activity of fission products with longer 
half.lives, although having lower fission yields, will begin 


to approach the activity of Zr-95, Nb-95, Pr-144 and Cs-134, 
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and the resolution of the system becomes considerably more im- 
pOoLctEanct® 

Zoeeioeene Lucure, the discriminator setting should be 
set higher (e.g. at 400 keV) than it was for these runs (52 
keV) to help reduce counting system dead times. Because it 
was desirable to minimize the need for movement of the detect- 
er, the first fuel element was placed in the holder prior to 
placing the dewar containing the detector on the top carriage. 
This made it impossible to make background runs prior to ob- 
taining the original spectra, and peaks noted in the channels 
below 900 were at one time Lisle: tomers 1Oiltacalt as mol DsSe — 
quent background runs proved this original assumption to be 
incorrect. With MITR fuel element 7M9 in the holder but with 
fe epeel 7mator, Cube memoved Gead times on the order of 15% 
were noted for a discriminator setting of 1.7, corresponding 
Mmemmemeui: OLL energy of 52 keV. For a discriminator setting 
of 4,2, corresponding to a cut off energy of 350 keV, the dead 


times were reduced to 33%, 


3) With the reduced background that will result from 
the less active spent power fuel pins, it would be possible to 
use a larger detector, which with its improved intrinsic 
efficiency will increase the number of counts measured, 
especially for the higher energy gamma rays. The use of a 
lacie nmEceteerortechGuld also help reduce the Compton edges 
somewhat by enabling more gamma rays to deposit their full 


energy. 
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4) Calibration of the system for the double escape peak 
of the 2186 kev Pr-114 gamma ray should also improve the deter- 
Mination of the Pr-144 content of the spent fuel. However, it 
MiOulLduee Motea the use Of a larger detector may decrease the 
probability of both gamma rays escaping from the detector after 
a Pair production event and thus decrease the size of the 
Seuple e€scape peak. 

In short, it is concluded that the proposed power fuel 
experiments using the apparatus described in this chapter are 
quite feasible. Sufficient Seal strength is easily obtain- 
Pemmonly 716th of the full collimator bore was used in the 
present scoping studies. Room background can be reduced to 
acceptable levels by proper choice of detector See ad 
Gebeeretiaekness and Shielding. Dead time can be reduced by 
musing a sufficiently high low level discriminator cut off. 
five kKeyelimitatiaon 1S really in the Compton background pro- 
duced in the detector. 

So far, then, we have investigated and established the. 
Beecential cbility to handle spent fuel and carry out a spec- 
EEOMeter Scan With available apparatus. In the next chapter 
a projection will be made concerning the suitability of the 


Gata obtainable with this apparatus for the fuel in question. 
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IV: CALCULATION OF GAMMA RAY INTENSITIES EXPECTED IN 
SEEN TEU 


fe introduction 

In order to determine the feasibility of studying the 
spent fuel made available by the AEC it was necessary to 
determine the expected intensities of those gamma rays use- 
ee LOT analysis. It was also necessary to determine whether 
they would be resolvable for a cooling time of approximately 
2.75 years. In addition to calculating the gamma rays emitted 
by fission products of U-235 fissions, an estimate was made of 
the expected gamma rays from fission products of Pu-239. 


weemecalculation of Expected U-235 Fission Product Gamma 


Ray Intensities. 


Calculation ae the expected intensities for gamma rays 
from U-235 fission products was performed by utilizing the 
calculated peak area (expressed as counts) obtained in Chapter 
3. The data obtained for MITR fuel element 5M27 was extra- 
._polated, using the familiar experimental decay Bee once 
to a 2.75 year cooling period, which is the expected cooling 
period for the Dresden fuel pins. The calculated intensities 
are shown in Table 13. 

Iimonic@eatemchieck tne Validity of extrapolating the 
intensity of a gamma ray, aS measured by the area under its 
peak, on the gamma ray Spectra, a parallel calculation was 
performed extrapolating the measured intensity for the gamma 


rays of 7M9 (cooling time = 2.75 months) to a cooling period of 9 


Sas 





Pepe tS: INTENSITY OF U=235 FISSION PRODUCT GAMMA RAYS 
AFTER 2.75 YEARS COOLING 


FISSION ENERGY PEAK 
PRODUCT (kev) AREA 
Pr-144 2186 22058 
1490 174.5 
1164 992.4 
696 G22 
Cs-134 796 A ones 
605 4,327.4 
Nb-95 . 766 93.3 
Zr-95 758 32.9 
724 19.9 
Cs-137 | 662 29,565.7 
Rh-106 623 544.0 
511 1280.0 


AG 





months. Corrections were made for the differences in counting 
prme, semickness Of tead sheets in the beam, and burnup. These 
calculated intensities are compared with the measured inten- 
Sity (peak area) for the gamma rays of 5M27 (cooling time = 9 
months) in Table 14. 

The calculated intensities compare well with the measur- 
ed intensities with the exception of Cs-137. No reason for 
the discrepancy could be discovered since both 5M27 and 7M9 
had similar irradiation histories and were used in Similar 
positions in the MITR core. 

In order to determine if the values of the calculated 
intensities listed in Table 13 were sufficiently strong to be 
resolved, their relative intensities were analyzed. This 
method of analysis was chosen because with spent fuel gamma 
spectra the difficulty is not one of obtaining sufficient 
count rate but rather of resolving lower energy gamma rays in 
the presence of Compton background from higher energy gamma 
rays. Since the height of the Compton edge and continuum are 
proportioned to the size of the peak it was felt that the 
ratio of a gamma ray's peak area to that of the sum of the 
peak areas for higher energy gamma rays would yield sufficient 
information to predict which gamma rays would be resolvable. 
Uiemresults Of tails analysis are Shown in Table 15. For com- 
parison the results of the GAMANL (16) calculation for fuel 
element 5M27 were analyzed in this manner and the results are 
leoteauin Column 3 Of Table 15. From this it was concluded 
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TABLE 14: CALCULATED AND MEASURED INTENSITY OF U-235 FISSION 
PRODUCT GAMMA RAYS AFTER 9 MONTHS COOLING, 


CALCULATED 
PEAK. 
AREA 


1,836.9 
LOS 43 
6,904.4 


8,880.6 
178,525.2 
69,640.8 
46 ,819.0 
6,980.4 
53,573.6 





Aa 


MEASURED 
PEAK 
AREA 


1914.0 
1034.5 
5884,4 
8003.9 
228,593.6 
80,577.3 
48 ,988.5 
6,535.3 
30,966.5 





Mee Ges 15; CALCULATED VALUES OF RELATIVE INTENSITY FOR U-235 
Piloo len me RODUCES APITIR 2,75 YEARS COOLING, 


CALCULATED MEASURED 
RELATIVE RELATIVE 
SiN OE Neier, TNTENS ITY 

(9 mos. cooling) 
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that all gamma rays were capable of being resolved by use of 
GAMANL (16). However, it is felt that the low absolute in- 
tensity of the peaks at 766, 758 and 724 keV would not lead to 
reproducible results in the calculation of the areas under 


the peaks. 
4.3 Calculation of Expected Pu-239 Fission Product 
4.3.1 Gamma Ray Intensities 


In order to obtain an estimate of the number of 
Pu-239 fissions tens have occurred in the fuel pins an energy 
balance was performed using power history data supplied by 
General miseeris for the fuel in question. It was assumed 


that the energy was released either by thermal fission of U- 


Meoror Pu-239, 1.é€.: 
E(TOTAL) = E(Pu-239) + E(U-235) 


Assuming 200 Mev are released per fission: 


Ze ae 
15,3g U-235 x Orzo x atoms , _gm mow > - l £ission 
. gm - mole 230 eGm U—2 35 1 atom 

: 200 Mev - Saige fission U-235 . watt - sec 


fission 683D fission + capture U-235 6.25 x 1014 Mev 


day Mid 15.3 (6.023) (577) (200) _ 45 95 my 


8.64 x 104 sec 10° watt 239 (683) (8.64) (6.25) 


Thus: 


E(Pu-239) = E(TOTAL) - E(U-235) 
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eC oes 2 20 


= 5,08 MWD 
TRO eA LOCOE fissions in Pu-239 fs 5 OS ite 
fissions in U-235 ee 


Since approximately 42% of the fissions are in Pu-239 
there iS a possibility of resolving Pu-239 fission product 
gamma rayS in the spectrum, The value of 42% is in good 
agreement with the values predicted using cross section data 
aS reported in Appendix A. 

From Table 8 it 1s seen that the fission yield of Rh-106 
- Ru-106 is increased from 0.4% for thermal fission of U-235 
to 4.4% for thermal fission of Pu-239, Table 13 shows the 
predicted peak area for the 623 keV peak of Rh-106 is 544.0 
counts for fission of U-235, The peak area for fission of 
Pu-239 would be 


4, 
4 


i 
oe 





> x BA 4440 (0 recounts 


oo 


The superimposed peak area for the fuel pins for the 623 keV 
i 
gamma ray would then be 
~416(5440.0) + .584(544.0) = 2574.0 counts 


This is an appreciable percentage of the peak area for Cs-137 


and thus should be resolvable, 
Ae Conc luslons 


A sufficient number of gamma ray peaks should be resolv- 


ol soiemotdtnotrecaliy significant 1 the spectra to enable 
Gs ileus 





nondestructive testing to be performed, Only one gamma ray 
need be resolved for each fission product. Even though the 
gamma rays of Nb-95 and Zr-95 will not be resolvable, the 
gamma rays of Rh-106 at 623 keV should be resolvable and thus 
the method of analysis utilizing three fission product activ- 


ity ratios recommended by Sovka (5) should be possible. 


me 





VC ONCLUSIONS 


The major conclusions reached as a result of the work re- 


ported in this thesis are as follows: 


iT) 


2) 


3) 


4) 


The calculated dose rate for the Dresden power re- 
actor fuel pins is well below that presently experi- 
enced from spent MITR fuel elements, and thus pre- 
sents no special handling problems. 

Sufficient spectral information should be obtained 

from the spent Dresden fuel pins to enable non-destruc- 
tive analysis to be performed. In fact, the increased 
energy resolution offered by present day Ge(Li) detect- 
ors and the 4096 channel analyzer coupled with the im- 
proved analytic capability of the GAMANL code (16) 
allow improved determination of the actual gamma ray 
mntensities Lrom Spent fuel spectra as compared to 
previous studies on MITR fuel by Sovka (5). 

Although Nb-95 and 4r-95 gamma rays will probably not 
be feeoived in the subject Dresden fuel, the expected 
relative iOEEae in Ru-106 and Rh-106 will permit sub- 
SeleutLon Of this, latter palr in the analytic procedure. 
A method of Compton suppression would be extremely 
valuable because of the intense Compton background pre- 
sent between 600 and 800 keV, the energy range for all 
but one of the gamma rays required for the method of 


analysis presented by Sovka (5). In addition, Compton 


23s 





4), “cone, 
Suppression might enable weaker intensity gamma 
rays potentially useful for analysis to be re- 
solved. The methods utilizing anti-coincidence 
would be extremely difficult to apply because 
of volume and weight limitations imposed by the 
scanning equipment. A duode type Ge(Li) detect- 
or (17,18), however, would be ideally suited to 
scanning of spent fuel Since it requires no 
additional space or volume than that needed for 
a standard Ge(Li) detector. 
As a result of this work, it can be said that all necessary 
groundwork for spent power reactor fuel pin analysis has been 
completed to the point where shipping and receipt of the fuel, 


and initiation of scanning experiments is now warranted. 


5 





APPENDIX A 


Pepaeamici=ealeulatvon of the ratio of the number of fission 
in Pu-239 to the number of fissions in U-235 was carried out us- 
ing one group theory utilizing effective neutron cross sections 
to account for the fission caused by epithermal neutrons and to 
mecount for all neutron captures in U-238. 

In order to obtain an expression for the number of Pu-239 
fissions we must first obtain an expression for the number of 
eee atoms converted to Pu-239 atoms as a function of time. 


For the one group model we can write: 


28 


dN 28 28 
dt aT On y ON (1) 
where: 28 ; 
N = no. of U-238 atoms at time t 
ey = efféctive microscopic capture cross section 
f 
EOr (U2 38 


® = one group neutron flux 


Solving Eq. (1) 


ZOO 28 
N = No lew-qy ey (2) 
where: 28 
No = initial no. of U-238 atoms 


For Pu-239 atoms we can write 


49 
an“? eee. 8 494.49 
= ony ¢X - GP so 
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N = no. of Pu-239 atoms at time t 


effective microscopic cross section for Pu- 
239 


» | 


Bestituting Eq. (2) for ee aaah) 100; Ae IS, 


49 


aN - 28 28 Ze 49 _49 
Ct — On, b Ns eo, pe ~V, . a 


Solving Eq. (4) yields 


28 
i = if exp ap olen Lo exp (-G'g t)- 


Ys 49 
Oe 
49 
exp(-(_ “@t) 
. 49 49 28 28 
Since N Oe ene and exp(- ¢_ E) 2 1 for norma! 
O a n,Y n,y 


fuel burnup the above expression can be approximated by 


28 
ae Gay eae 7-"*¢ t) coy) 
© 49 a 
O, 


For the time ratio of Pu-239 fissions we can write 


49 
aqF-  -491 49 
a OF WF DN (6) 


peDstaturcinGg Eq. (5) for see in Eq. (6) we have 


eG 
ar" = Op PN | (n,y eee en( (-°7¢ t) 
G49 
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Integration yields 


28 
49, 28] On 49 49 
Fa = CF pN inti || = l-exp (¢ D t) 7 p (7) 
O a a 
Gg 49 
a 
( 
Expanding the exponential term of Eq. (7) in a Taylor Series 
yields 
i 28 
49 _ 49 28 inl 59% 49 _ 49 2 
F = ¢ PNS a a ee ote /2 ti) rely 


Ignoring all terms above the ae term yields 


28 
coe G96 n?8 (ny | ero es (9) 
6-49 : 
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For the number of U-235 atoms we may write 


25 | 
aN a 2g aS 
; no N (9) 
where: 25 
N =no., of U-235 atoms at time t 
25 : . 
Om = effective microscopic cross section for 
U-235 
Solving fq. (9) 
25 25 Zo 
N aes exp (- QD t) eno) 


Sea 





' The rate of U-235 fissions iS given by 


ZS 
aF 25 Bs 
ae = OPN (11) 
where: 25 
F —meEoeatlemo., of U=-235 £LlSsions £from time £ =0 
Pmoetieutcing Eq, {10) into Eq. (11) 
ar 7 _ 25 25 
7 Buo*? (1-exp(- GP t))/G& (12) 


The ratio of Pu-239 fissions to U-235 fissions using Eq. (8) 


and Eq. (12) becomes 


49 28 28 49 2 25 
e9 ge \ING ony (A/G PEGS D 
eee ioe «seme tayo 
F G49 ace 6225 1- ay? SM, 25) 
a iE 
49 28 28 
- Gg No ny 23 
- Fomor (pt) (13) 
- eS 25 
is ae hy (5 


In order to eliminate flux and time dependence of Eq. (13) 


return to Eq. (10) and take logarithms of both sides 
In (N* > no? = ee 
Solving for gt yields 
gt = a7G >) An 29 ye? 


om 





PUSSELEMEING LOm et in Eq. (8) yields 











25 2 

49 28 Ze N : 
mee ay Ge Ga Ne ibe 2 C14) 
i aa ee 25 

Zo 25 ZS 2° 25 N 
F Og (5 Ne ~ N 
Rearranging, yields 

49\/ 28 28 2G ele 
po? £ ba . lt we 
——— = 1/2 | ——— 2 STE ee era 
pe? g°° ia Z5 1- (N?? yn 7° NZ? 


The ratio on, , Gn? is defined as the initial conver- 


Sion ratio (ICR) and is equal to the initial number of neutron 
absorptions in U-238 divided by the number of neutron absorp- 


tions in U-235., Thus the final working equation becomes 





2 
49 25 
a ae ; = ) ICR ve (15) 
F e (a-n?5 ny 25) es 


In order to apply Eq. (15) we need to Specify the erroce 


tive cross sections via 


mal and epithermal events the following approximate formula- 


and GC : To account for both ther-— 


tion can be employed: 


C = bh + Si caylee |* Re 
£ 


where: J = effective microscopic cross section 
2510s 





Gen = ene rina: eee Duc Gross section 
G55 = measured ratio of epithermal to thermal fissions 
RI = infinite dilution resonance integral 


Using tabulated values for G49) and RI(20) and a value of 
e Ze eee Piecaliavalue £or the lattice in question, we 


obtain GE = 826 barns and ee = 660 barns. Using these 


ol 
values for G and an omy Omecdime verkiye ror ECR Of 05, 
and the values of N* and te? from Table 1 we obtain from 
Eq. (14) a value of ae = O02 


Ge 
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